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Abstract. We present a uniform analysis of the ROSAT HRI and the ASCA GIS/SIS data for 79 distant clusters of galaxies in 
the redshift range 0. 1 < z < 0.82 to study the global structures of the intracluster medium. We have constructed an X-ray catalog 
consisting of the largest sample of clusters in the redshift range for which pointed X-ray observations were carried out with 
both the observatories. We determined the emission-weighted X-ray temperatures of the clusters with ASCA, while we studied 
surface brightness distribution with the ROSAT HRI utilizing the isothermal /? model. We investigated the statistical properties 
and trends for redshift evolution of the X-ray parameters including the temperature, the density profile of the intracluster gas 
and the gas-mass fraction within rsoo- We also present correlations of the cluster parameters with the X-ray temperature and 
with the core radius and compare them with the predictions of the self-similar model, from which we discuss the possible origin 
of the double structure discovered in the core radius distribution. 
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1. Introduction 

Clusters of galaxies are the largest collapsed systems known in 
the universe. Because the time scale of evolution of clusters is a 
significant fraction of the age of the Universe, the clusters may 
preserve information about the early universe and thus are con- 
sidered to be excellent tracers of the formation and evolution of 
structures. They are considered to continue to grow into larger 
systems through complex interaction between smaller systems, 
namely merging process. It is possible that the clusters that we 
see are in different stages of evolution. 

For the understanding of cluster structure and evolution, we 
believe it is important to analyze systematically a large number 
of cl usters at various re dshifts. After the ASCA jTanaka et all 
1 1994 and the ROSAT jTriimpeit Il993l) X-ray observatories 
were put into orbit, it became possible to study relatively dis- 
tant clusters at X-ray energies. During 7-10 years of observa- 
tions, more than one hundred clusters were recorded with both 
observatories. Since ASCA has a high sensitivity to measure 
the X-ray spectrum in the wide energy band while ROSAT is 
good at imaging in the soft X-ray band, the two observatories 
are an excellent combination to study properties of the intra- 
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cluster medium (ICM). At present the XMM-Newton and the 
Chandra satellites are in orbit and generate much cluster data 
with higher sensitivities. However, as we will mention below, 
the data set used in the present paper will be one of the best 
existing to construct the largest sample of distant clusters and 
study global X-ray structures. 

X-ray observations bring us valuable information on not 
only the physical state of the ICM but also the underly- 
ing potential structure of the clusters. Statistical studies are 
very powerful in exposing the physical nature of the clus- 
ters. In particular, nearby clusters have been extensively stud- 
ied a t X-ray wavelengths. As to t he low redshift samples (z < 
0.1). iMohr. Mathiesen & Evrardl Jl999) performed a system- 
atic analysis on the ROSAT PSPC data of 45 clusters and utiliz- 
ing the published ICM temperatures, they investigated the cor- 
relation between the ICM mass and the temperature, namely 
the M gas - T relation. They found that the slope of the relation 
is significantly steeper than that predicted from the self-similar 
model jKaisen Il986l) . The observed luminosity-temperature 
(Lx - T) relation is also known to be steeper than the e xpec- 
tation of the self-similar model (e.g. Davi d et all Il993l) . The 
inconsistency between the observations and the simple theoret- 
ical model has been debated for many years and various possi- 
bilities such as non-gravitational heatin g (e.g. E vrard & Henrvl 
1991t lCavaliere. Menci. & T ozzi. 1997) and dependence of gas 
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mass or gas-mass fraction on the temperature have been pro- 
posed (e.g. lDavid et al.lll993llNeumann & Arnaudll200ll) . For 
higher redshift samples, however, the number of clusters that 
were uniformly analyzed was limited compared to the nearby 
clusters. The ASC A spectroscopic data of d istant clusters were 
compiled by Mus hotzkv. & ScharJ i 19971) (38 clusters with 
z > 0.141 iWhitd feOOOl) ( 41 clusters wi t h z > 0.1 and 65 
clusters with z < 0.1), and iNovicki et al 1 d2002l) (32 clusters 
with 0.3 < z < 0.6 and 53 clusters with z < 0.3), while the 
combined (i.e. spectral and imaging) analyses of the distant 
clusters were separately published in lAllenl i 19981) (13 clus- 
ters with 0.1 < z < 0. 451 iHashimotodanild 19991) (27 clusters 
with 0.1 < z < 0.781 lEttori&Fabianld 19991) (3 6 clusters at 
7 >_0.05l lSchindlerl ll999) ( 1 1 clusters at z > .3llLewis et al 



( 1999) (14 clusters with 0.14 < z < 0.55) and Vikhlinin et al 
1 20021) (22 clusters at z > 0.4). Recently, lEttori et al. I d2004 
reported the Chandra analysis of 28 clusters at 0.4 < z < 1.3 
and observed the steeper slopes in the L-T and the M gas - T 
relations, which provided hints of negative evolution in their 

relations at high redshift 

On the other hand, |Qti J200ll) and lota & Mitsudal d2002l) 
have attempted to construct the largest X-ray sample of dis- 
tant clusters with z > 0.1 based on the combined analysis of 
the ROSAT HRI and the ASCA GIS/SIS data. The data set used 
in the analysis is one of the best suited for the purpose of the 
present systematic study for the following reasons: the largest 
cluster sample of the pointed observations were stored in the 
ROSAT and the ASCA archival database, which enable us to 
cover the widest redshift range up to z ~ 1 . The focal plane 
instruments, the ROSAT HRI and the ASCA GIS/SIS have suf- 
ficient sensitivities to study the global properties of the ICM 
spatial structure with a typical resolution of 5" and spectral 
features such as the ICM temperature and the luminosity, re- 
spectively. In addition, the instruments' responses were well 
calibrated and the use of the same detectors for all the clusters 
will reduce systematic effects when comparing their results. 
Our sample covers the intermediate redshift range of 0.1-0.8 
between those of Mohr et al.'s sample (z < 0.1) and Ettori et 
al.'s sample (0.4 < z < 1.3), thus combining our data with 
the other two samples is helpful in discussing the evolution of 
internal structure of the ICM directly from the observational 
point of views. 

In this paper we describe a uniform analysis of the ROSAT 
HRI and th e A SCA SIS/GIS data of 7 9 clusters performed in 
Otal J200ll) and lOta & Mitsudal J2002) and thus provide an X- 
ray database of ICM structure with the widest redshift range 
of 0.1 < z < 0.82. In order to perform a reliable parameter 
determination we paid special attention to evaluating all pos- 
sible systematic errors in the spatial and the spectral analyses. 
We determined the spectral and spatial properties of the intr- 
acluster gas in terms of the temperature, the core radius and 
the central electron density etc. for all of the clusters. In the 
image analysis we found that the sample can be naturally clas- 
sified into regular and irregular clusters according to the X-ray 
surface brightness distribution. We present the statistical prop- 
erties of the X-ray parameters and the relations to the double-/? 
nature of the clusters discovered in the core radius distribution. 
We also estimated the cluster total mass, the gas mass and the 



Fig. 1. (a) Redshift distribution of 79 distant clusters, (b) 
Redshift distribution of the 45 regular clusters (open) and the 
34 irregular clusters (hatched). 

gas-mass fraction within a radius corresponding to a fixed over- 
density of 500. We further studied the redshift dependency of 
the parameters and their correlations, where we considered a 
systematic error due mainly to the limited sensitivity of the in- 
struments and some assumptions used in the estimations. 

This paper is organized in the following manner. In the next 
section, we describe the characteristics of our samples. In sec- 
tion |3 and 0] we explain the methods of spatial analysis with 
ROSAT and spectral analysis with ASCA in detail. In section[5] 
we derive the X-ray parameters to describe the global structure 
of clusters and investigate their redshift dependence. In section 
[6]we study correlations of the parameters and discuss the prop- 
erties of the ICM. In sectionQwe summarize our results. 

We use Qm = 0.3, Qa = 0.7 and hjo 
Us, 770 km s^'Mpc -1 ) = 1. The quoted errors are the 90% con- 
fidence range throughout the paper except where noted. 

2. The sample 

We have selected distant clusters with 0.1 < z < 1 for 
which pointed X-ray observation data with both ASCA and 
ROSAT/HRl are available. Though there are 83 clusters which 
meet the criteria, we rejected three (A222, A223,A1758S) be- 
cause of the large (> 210") pointing off-axis angles in the 
ROSAT observations. We did not include a high-r edshift clus- 
ter, AXJ2016+112 at z = 1 faattoriet all 1 19971) because the 
Chandra observation showed that the X-ray emission from the 
direction of AXJ2016+1 12 can be resolved into point sources 
and th e cluster diffuse emission is not prominent dChartas et all 
120011) . The final sample comprises 79 clusters. Among them, 
three clusters (#1 PKS0745-19, #6 A2204, an d #13 A1689) 
overlap with lMohr. Mathiesen & EvrardMigggh 's nearby clus- 
ter sample and thir ty are known to hav e strong gravitational 
lensing sig nals (e.g. lHattorie7ail ll999). The observation logs 
are summarized in Table 1 , where the clusters were sorted ac- 
cording to the redshift and numbered from 1 to 79. The redshift 
distribution of the sample is shown in Fig.[^. The clusters with 
0.1 < z < 0.3, 0.3 < z < 0.5 and 0.5 < z < 1 make up 58%, 
32%, and 10% of the sample, respectively. The average redshift 
is 0.30. 

Because our analysis targets were collected from the pro- 
posal observations and the sensitivities for high-redshift clus- 
ters are limited, we have to carefully consider possible selec- 
tion bias. In the first step, we compare the sample with other 
unbiased cluster samples. We show temperature distributions of 
our sample and t he nea rby X-ray flux-limited 55 cluster sample 
constructed by Ed ge et alJ dl990) in Fig. [2] Our sample covers 
the equivalent temperature range, but has a higher average tem- 
perature of 6.8 keV. The Kolmogorov-Smirnov (K-S) test gave 
the probability that the two samples are from the same tem- 
perature distribution as 0.06 (the K-S parameter, D = 0.24). 
Observation bias will be discussed in fl5.1.1l and l5.2.2l in more 
detail. 
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Fig. 2. Comparison of the temperature distribution of our 79 
dista nt clusters (o pen) to the nearby flux limited 55 sample of 
iFxtee etal J ill 99(1) (hatched). 

3. Spatial analysis 

3.1. Data reduction from event lists 

We have retrieved the event lists of the ROSAT HRI detec- 
tor from the ROSAT Data Archive of the Max-Plank-Institut 
fur extraterrestrische Physik (MPE) at Garching , Germany. We 
used the EXSAS analysis package (Zimmerma nn et all 1 1992) 
to produce X-ray images from the event lists. The raw HRI data 
has a 0".5 spatial resolution. However, since the half power di- 
ameter of the X-ray telescope is 4". 8 at the optical axis and 
the photon counting statistics of the present data are limited, 
it is not worth oversampling the telescope point spread func- 
tion (PSF). We thus rebinned the image into 5" bins, within 
whose diameter about 70% of photons from a point source are 
included. 

If there were multiple pointings for a cluster, and the point- 
ing offset angles between them are smaller than 3', we super- 
pose the event lists in the sky coordinate. The PHA channel 
of 1-10, which corresponds to 0.2-2 keV, was used to avoid 
particle background events. 

Next we searched for foreground/background sources in the 
field of view with the standard source detection program in the 
EXSAS. We created lists of all the sources detected by the max- 
imum likelihood technique and having a likelihood of > 8. We 
used the lists to exclude those sources from the analysis. 

3.2. Centroid determination and X-ray morphology 

We apply a method to determine cluster centroids and clas- 
sify the morphology of clusters. A similar techn ique to fin d the 
cluster center was first introduced by Mohr, Fabricant & Geller 
dl993h . We extended the method to evaluate the systematic er- 
rors of the centroid determination for distant clusters and found 
that the analysis can also be used to determine the regular- 
ity of the X-ray distribution. The cluster emission is typically 
extended about 4' from the optical axis of the telescope. In 
this image region, the position depend ence of the telesco pe vi- 
gnetting is known to be less than 2% dBriel et all Il997^. The 
position dependence of the background intensity is also small 
up to ~ 12'. Thus they do not affect the centroid determination. 

The analysis consists of two major steps. First we estimate 
the 0-th order center position and a parameter which represents 
the extent of the X-ray image, utilizing 1 -dimensional image 
projections and Gaussian fits to them. Then, starting from the 
0-th order initial value, we determine the center from the center 
of gravity of the photon distribution. 

First, we extracted a 4' x 4' image that contains the cen- 
tral region of the cluster emission. Then we projected along 
the x- and y- axes to get one-dimensional intensity profiles. 
We fit these with Gaussian functions, and determined the 
mean, (xc,o,yc,o) an d the width, (cr A . o, <x v o). To determine these 
values with higher accuracy, we extracted an image of size 



3o~ x$ x 3cr v o whose center is at (xc,o,yG,o) an d performed 
the one-dimensional Gaussian fitting again to derive the next 
set of (xc,i,yc,i) and (cr A% i, cr v i). The procedure was iterated 
; times until the mean position converged within 0.1 pixels 
i.e. \(x CJ ,y CJ ) - (x G j-i,y GJ -i)\ < (0.1,0.1) or the number of 
iterations reached i - 20. We then define a parameter <x as 
o~ = (<t x j + <x v ,)/2, which is a measure of the image extent for 
further analysis. 

Next we determine the cluster center from the center of 
gravity of the photon distribution in an aperture circle of radius, 
R. Then if the X-ray image is circularly symmetric and the cen- 
ter of the circle is at the X-ray center, the center of gravity of 
the photon distribution should coincide with the center of the 
circle. For a given value of R, we can determine the center po- 
sition, r,- = YirPI Er 1' where p represents the position of the 
photon, and 2r sums all the photons over the circular area of 
the radius R. Then starting with the mean position determined 
in the previous paragraph, tq = (xcj,ycj), we extract a circular 
image of the radius R centered at ro, and calculate the centroid 
position r\. We continue the iteration until |r,- - r,_i| becomes 
less than 0.1 pixels. If there were contaminating sources in the 
circle, we excluded the region centered at the sources and the 
region symmetric to them so as not to affect the centroid deter- 
mination. We changed the aperture radius R from 2 <x to 9 & to 
study the R dependence of the results. As a result, while some 
clusters showed constant centroids almost independent of the 
radius, others exhibited systematic behavior. Two representa- 
tive cases are shown in Fig. [5] 

The dependence of the center of gravity on the aperture 
radius may indicate some asymmetry or substructures of the 
cluster. However, because the results with different R are not 
statistically independent, the Poisson fluctuations may produce 
systematic deviations. To estimate this effect, we performed 
Monte-Carlo simulations and generated a number of simula- 
tion images assuming the isothermal jfi-model as the X-ray im- 
age distribution. We calculated the RMS value of the distance 
between the center of gravity and the true center, cr M c, for var- 
ious combinations of parameters, the cluster core radius r c , the 
number of photons N and the size of the aperture R, while f3 
was fixed to the typical value of 0.67. 

If the displacement of the center of gravity with different 
values of the aperture radius R is larger than the standard de- 
viation determined from the simulation, we can conclude that 
there is asymmetry or substructures. In Fig. [5] we show the de- 
viation of the center of gravity from that determined for R = 3dr 
normalized by the standard deviation. For some clusters the de- 
viation is well within the lcr M c level; however, for others it is 
not. We define criteria for irregularity of clusters as follows: (1) 
more than 1 data point whose deviation exceeds 3cr M c, and/or 
(2) deviations (> lcr MC ) seen in more than a few consecutive 
points. Accordingly, we classified our sample into regular and 
irregular clusters. The results are listed in Table 2. The ratio of 
regular to irregular is 45 : 34 and the redshift distributions of 
the two subgroups are shown in Fig.^i. 
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Fig. 3. Centroid determination and classification of X-ray morphology. #13 A1689 and #19 A2163 are shown in the left and right 
panels, respectively. From top to bottom: (a) the HRI image in the pixel coordinates, (b) The centroids determined in circles of 
radii R = ncr (n = 2, 3, .., 9), where the size of the symbols are nearly proportional to n. (c) The deviations the centroids relative 
to that determined for 3<r in units of standard error, shown as a function of R [&]. While A 1689 is classified as a regular cluster, 
A2163 is an irregular cluster due to the significant centroid deviation. 



Fig. 4. Radial X-ray surface brightness profile of #13 A1689 fit- 
ted with the single y6-model (a) and the double /3-model (b). The 
crosses denote the observed radial profile of the ROSAT HRI, 
and the step functions show the best-fit y6-models. The best-fit 
background levels are shown with the horizontal dashed lines. 
In the panel b), the contribution of the inner and the outer com- 
ponent of the double /?-model are also shown with the dash-dot 
and the dotted lines, respectively. 

3.3. Radial X-ray surface brightness profiles 

We derive azimuthally-averaged radial profiles of the X-ray 
surface brightness centered at the cluster centroids that were 
determined within the 3cr aperture radii for both the regular and 
irregular clusters. We have chosen the centroids because they 
are less affected by the Poisson fluctuations in the outer image 
regions. The bin size of the radial profile is 5". We excluded the 
contaminating sources from the integration area with a circle of 
radius 5 times the FWHM of the PSF at position. 

We study the effect of the choice of the center on the re- 
sults of the radial profile fits, by shifting the center positions 
determined for R between 2cr and 9<x. We found the yS-model 
parameters ( 33. 4t constant within the statistical errors for all 
the regular clusters and many of the irregular clusters. For 30% 
of the irregular clusters, the resulting core radius increases typ- 
ically by a factor of 2 as R increases from 2<x to 9<x. 

3.4. Radial profile fitting with a single fi-model 

In the following two subsections, we analyze the 
X-ray radial profile und e r the isothermal yS-model 
JCavaliere & Fusco-Femianol Il976h . The single yS-model 
fitting function is written as 



S(r) = S c 



1+1^ 



-3/3+1/2 



+ c, 



(1) 



where So, r c and j3 are the central surface brightness, core ra- 
dius and the outer slope, respectively, and C is a constant back- 
ground. In Fig. |4k we show an example of the radial profile 
fit. 

It is crucial for reliable determination of the model parame- 
ters to estimate the background level correctly. Th e HRI back- 
ground is dominated by the particle background (Briel et al., 
1 19971) and the detailed calibration by Snowden ( 1998) showed 
that it can be modeled by a constant image within ~ 12' from 
the detector center. The counting rate of the particle back- 
ground depends strongly on the satellite orbit and time, and 
typically ranges from 1 to 10 counts s over the entire detector. 
Thus the background counting rate needs to be determined ob- 
servation by observation. We determined the background level 



Fig. 5.x 2 contour of the single /3-model fit for #13 A1689. The 
position ofx 2 minimum is denoted with a cross. The curve cor- 
responds to the single-parameter error domain at 90 % confi- 
dence. 

Fig. 6. Effect of the outer cut off radius in the single /3 model fit- 
ting. The results for #13 A1689 are shown as an example. The 
four parameters of the single y6-model, S o [counts s~ 1 arcmin" 2 ] , 
P, r c [arcsec], C [10~ 3 counts s _1 arciriirr 2 ] are shown as a func- 
tion of the outer cutoff radius, r out in unit of cr. The results of 
the constant fits to the data points at r out > 10 are shown with 
the dotted lines. 

from the radial profile including the background as one of the 
fitting parameters as Equation^ Since we do not know the true 
extent of the cluster emission and also r c and /3 are strongly cou- 
pled (Fig. the background level and the model parameters 
depend on the outer radius of the fitting area, r out . Particularly 
when r out is too small, the background level determined from 
the fit tends to be over- or under-estimated and results in un- 
certain r c and (3. However, as shown in Fig. |6j the best-fit pa- 
rameters converge to constants if r out is large enough, typically 
> 10cr. For all the clusters, we confirmed that the parameters 
converge at r out = 12'. Thus we decided to adopt this value 
for all clusters, within which the background can be regarded 
as constant. Finally the background level we obtained from fits 
are reasonable background levels. 

Because the effective area of the X-ray telescope decreases 
with off-axis angle, the cluster image may be affected by the 
vignetting. However since the off-axis angle of the cluster cen- 
troid of the present data is smaller than 4', and the typical spa- 
tial extent of the present clusters is ~ 4', the vignetting effect is 
at most 3 % at the rim of the clus ters (the HRI vignetting func- 
tion is given in lBriel et all (fl997)'). Although we performed fits 
with the /?-model fitting function with and without the correc- 
tion of the vignetting function, the results showed no differ- 
ence. 

The /3-model function needs to be convolved with the X-ray 
telescope PSF, then integrated over the image bin. However, 
since the present image bin size is larger than the extent of the 
PSF, the convolution with the PSF is not important. Also, the f3- 
model function varies slowly within the 5" bins in most cases, 
so integration within the image bin can be replaced by the value 
at the center of the bin. We confirmed these with simulations 
by comparing two cases with and without convolution in the 
fitting model. As shown in Fig.0 we find the difference of the 
two cases to be negligibly small as long as the core radius is 
larger than the bin size, 5". We also confirmed that both cases 
well reproduce the assumed /3 value of 0.6 within reasonable 
statistical errors. Thus in order to make the computation time 
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Fig. 7. Reproducibility of r c in the single j6-model analysis. The 
jc-axis is the assumed core radius in the simulation cluster im- 
age, and the y-axis is the core radius derived from fitting with 
two different models: the yS-model with and without the PSF 
convolution. The results for three different photon counts are 
shown in each panel. 

shorter, we skip the convolution with the PSF and the integra- 
tion over the bin. If the best-fit value (and the error domain) of 
the core radius is smaller than the bin size, we regard it as an 
upper limit. 

In Fig. |5] we showed a^ 2 contour map on the ft - r c plane 
where the other free parameters, So and C, are optimized at 
each point of the plane. The two parameters are strongly cou- 
pled and the allowed parameter range exists in an elongated 
region. We quote the 90% confidence intervals on the best-fit 
parameters. Table 2 lists the results of single /i-model fitting. 
For #43 A1758 and #79 MS 1054.5-0321, because the fitting 
parameters did not converge properly, we optimized the param- 
eters within the range of f3 < 3. 

We compared th e results of /i-model fitting with 
iHashimotodanil ill 999t) for 27 RO SAT/UR1 clusters with 0.1 < 
z < 0.78 and lEttori et al.l J2004 for 10 high redshift Chandra 
clusters with constrained model parameters. We found that 
there is a good agreement for /3 and r c measurements between 
our and Hashimotodani samples within their statistical errors. 
Furthermore for 8 of the 10 high redshift clusters, there is 
a good agreement between our and the Ettori et al. samples 
within the 90% statistical errors. For the two highest redshift 
samples, #70 RXJ1347.5-1 145 (z = 0.451) and #71 3C295 
(z = 0.4641), r c is systematically smaller by about 40% com- 
pared to the Chandra measurements. /3 is also systematically 
smaller with the mean ratio of 0.93 for 10 high redshift sam- 
ples, although the difference is the ~ 2<x effect. We consider 
that the limited sensitivity of the HRI for the outer part of the 
cluster emission may cause the underestimation of /? for the 
highest redshift (z > 0.4) samples. 

As a result, the fractions of clusters with^- 2 values exceed- 
ing the 90% and 99% confidence levels are 25/79 (=0.32) and 
11/79 (=0.14), respectively. The fractions are larger than ex- 
pected only by the statistical errors. Thus there may be some 
systematic errors that are not well explained by the single /3- 
model. For some clusters, systematic deviations from the cur- 
rent model are particularly seen in the central regions (e.g. #13 
A 1689), which will be discussed in the next subsection. 

We also evaluated an X-ray significance radius, rx, repre- 
senting the extension of the observed cluster X-ray emission at 
which the best-fit y6-model surface brightness becomes equal to 
the 3cr background level. The results are also listed in Table 5. 

3.5. Radial profile fitting with a double fi-model 

For some clusters, systematic residuals are seen in the results 
of single yS-model fitting. As often seen in nearby clusters, 
this may be attributed to the presence of central excess emis- 
sion. The excess component is ofte n represented by an ad- 
ditional yS-model component (Jones & Format! Il984h . whose 



Fig. 8. Effect of the inner cutoff radius in the single /3 model fit- 
ting. The results for#13 A1689 are shown as an example. From 
top to bottom, the resultant values of r c [arcsec], F (see 33. 5 I f or 
definition) and the probability of exceeding the F value, P, are 
shown. 

core radius is 10 ~ 200 /ijjkpc, and on average 60 /i^jkpc 
(Mohr. Mathie sen & Evrardtll999l) . This corresponds to only 
several bins of the present analysis even at z — 0. 1 . It is not easy 
to constrain such central emission in the distant clusters. To 
evaluate the statistical significance of the central emission, we 
attempted two methods. We restrict this analysis to the 45 reg- 
ular clusters because the irregular clusters are often accompa- 
nied by substructures, which can cause artificial double struc- 
tures in their radial profiles. 

In the first analysis, we exclude central bins from the fit and 
investigate the variation of single /3-model parameters against 
the inner cutoff radius, n n . We test the improvement in the^ 2 
value of the fit with the F-test compared to the case of r m = 
(Fig- ID- We define the value of F as F = ((x\ - X^liyx - 
v 2))/(X2/v2), where x 1 and v are the minimum^ 2 value and the 
degrees of freedom, and the suffix "1" and "2" correspond to 
the case of r m = and r m = n pixels, respectively {v\ — V2 = «)■ 
For nine clusters, we find that the probability of exceeding the 
F value, P, rises at a certain r m . For seven such cases, the core 
radius also starts increasing at the same inner cutoff radius. This 
indicates that the nine clusters have significant two core sizes. 
Moreover for the seven of the nine the core radius of the inner 
component at which the significant improvement of x 1 starts is 
roughly r c . 

In the second analysis, we assumed the double yS-model 
composed of different core radii {r\ < r^), 



S(r) = ^Si 



1 + 



-3/3/+1/2 



+ c, 



(2) 



and performed fitting to the radial profiles (Fig. @J>). Because 
the inner slope j3\ is insensitive to the fit, we linked it to the 
outer value, fi\ = fii- We tested the statistical significance of the 
improvement of the \ 2 value by the F-test against the single f3- 
model; F = ((y 2 -^)/(v s - v d ))/(xl/v d ), where the suffixes "s" 
and "d" stand for the case of the single yS-model and the double 
yS-model, respectively. We found that nine of the regular clus- 
ters have significance above the 95 % level. Those nine clusters 
are the same clusters for which the double core nature is found 
in the previous analysis. We refer to the nine clusters with sig- 
nificant double structure as "double-/?" clusters and show the 
result of the double /3-model fitting in Table 3. 

We have plotted the two core radii against the core radius 
derived by the single /i-model for the nine "double-/?' clusters 
in Fig. |9^. It is remarkable that one of the two cores is nearly 
consistent with the core of the single /J-model, namely r c ~ r\ 
or r c ~ r2- This indicates that the single /3-model tends to de- 
tect the one of the two components that is more dominant. We 
also show the ratio of two cores versus the ratio of two normal- 
ization factors in Fig.|9j5. r2/r\ is ~ 4 on average, while S2/S 1 
ranges from 0.01 to 1 . These are consistent with nearby clusters 
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Fig. 9. Core radii derived by the double /?-model for the nine 
clusters with significant double structure. The panel a) shows 
the relation between two core radii of the double /(-model and 
that of the single /(-model, where the inner and outer cores of 
the double /? model are shown with the filled and open circles, 
respectively. The panel b) shows the ratio of the two cores and 
the ratio of the two normalization factors in Equation 13 The 
two data points whose S2/S1 are large (~ 1) in the panel b) 
correspond to the two clusters with large core radius of the sin- 
gle /(-model in the panel a). 

jMohr. Mathiesen & Evrardlll999l) . We can classify them into 
two groups: inner core dominated clusters (S2/S 1 < 0.1) and 
outer core dominated clusters (S2/S 1 ~ 1). For the inner core 
dominant clusters, the single /?-model fit picks up the inner /?- 
model component, while for the outer core dominant clusters, 
the outer component is picked up. In Table 2 the inner/outer 
core dominant clusters are denoted with W(l)/W(2). 

We then investigated the reproducibility of the best-fit pa- 
rameters of the double /(-model utilizing Monte-Carlo simula- 
tions. We made ten simulation clusters for each set of model 
parameters described below and fitted them with the double /?- 
model. We assumed r\ — 3 pixels(= 15") as a typical core 
radius of the inner component and several different values of 
r2 between 6 and 24 pixels. We fixed fi\ — ($2 — 0.67. For 
the intensity ratio, we checked two typical cases, i.e. S%/S 1 = 
0.1 and 1.0. We then generated 75000 events within a radius 
of 12', which are typical total counts for double-/? clusters in- 
cluding background. We confirmed that the two core radii are 
well determined as long as r%lr\ > 3. Thus, at least for the nine 
double-/? clusters we detected, we can conclude that the model 
parameters are well-determined by the model fitting. 

3.6. Fraction ofdouble-fi clusters 

We also find that the double-/? clusters are found only at z < 0.3 
and the ratio to the regular clusters at z < 0.3 is 32%. However 
we have to take into account the fact that such double struc- 
tures are difficult to find at higher (z > 0.3) redshifts because 
of the low surface brightness. In order to constrain the double 
structure at higher redshifts, we estimated the upper limits of 
the possible additional /(-model component. For this purpose 
we fitted the radial profile with double /(-model with the ratio 
of the two cores fixed at the average of nine double-/? clusters, 
ri/ri = 4. We need to consider two cases: (1) the inner compo- 
nent is dominant, namely the inner core radius is approximately 
the core radius of the single /?-model, r\ ~ r c , and (2) the outer 
component is dominant, r2 ~ r c . We thus performed the fit with 
r\ ~ r c or r2 ~ r c as the initial value. For some cases, the fit 
converged to certain best-fit values or upper limits of the fitting 
parameters. However in some cases the fit did not converge. In 
such cases we fixed r t or r2 at the value of r c obtained from the 
single /(-model fitting. 

We find that in addition to the nine double-/? clusters, there 
are several other clusters that may contain a second compo- 
nent. In case (1) the outer component is marginally detected 
for seven clusters at the 90% confidence level. For other 17 



clusters, the upper limit of the surface brightness of the outer 
component is within the range (0.01 - 0.1)5 1 which is compa- 
rable to the range for the seven inner-core dominant double-/? 
clusters. In case (2) the inner component is marginally detected 
for three clusters, and the upper limit is consistent with the two 
double-/? clusters for the other 16. In Table 2 the clusters with 
marginal inner/outer component are denoted with S(l)/S(2). 

We finally obtain the fraction of clusters with marginal 
double-/? structure at z > 0.3 to be 35%. It is comparable to 
the occurrence of the double structures in the lower redshift 
systems within the Poisson errors. Thus we cannot conclude 
that there is significant evolution of the fraction of the double- 
ts clusters in the regular c lusters within the o bserved redshift 
range. On the other hand, lEttori et al. I (120041) noted based on 
the /?-model analysis that their high redshift (z > 0.4) sam- 
ples do not show any significant double structure in the surface 
brightness distribution. Thus there might be a trend of evolution 
in the core structure of the ICM distribution starting around 
z ~ 0.4. We suggest that this should be confirmed by further 
observations. 

4. Spectral analysis 

4.1. Data reduction 

We retrieved the ASCA data sets from the High Energy 
Astrophysics Science Archive Research Center (HEASARC) 
at NASA/Goddard Space Flight Center and the DARTS Online 
Service at the Institute of Space and Astronautical Science 
(ISAS) that were screened with the standard REV-2 process- 
ing. We use the FTOOLS analysis package to reduce the cluster 
spectra, and calculate the response functions of the telescopes 
and the detectors. 

The GIS was operated in the PH-nominal mode during ob- 
servations. The SIS has several choices between the CCD data 
modes (FAINT or BRIGHT) and the CCD clocking modes 
(1CCD or 2CCD or 4CCD). For observations done in a mix- 
ture of the FAINT and BRIGHT modes, we can combine con- 
verted BRIGHT (on-board FAINT) and on-board BRIGHT 
mode data. 

First we extacted the X-ray images in the 0.7 - 10 keV for 
the GIS and the 0.5 - 10 keV for the SIS, respectively. The two 
sensors in the same system, namely GIS-2, -3 and SIS-0, -1 
are added together. We fitted the projected images to Gaussian 
functions and determined the peak positions. We accumulate 
spectra from a circular region centered on the Gaussian peak, 
where the extraction radii are 6' and 3' for the GIS and the SIS, 
respectively. We will discuss the effect of the contamination 
from foreground/background point sources in the next subsec- 
tion. We select a larger integration area for the GIS because the 
FWHM of the point spread function of the GIS detector alone 
depends on the incident pho ton energy E keV, and is giv en by 
0.5 x (5.9/F) dThe ASCA Data Reduction Guide IkOoA The 
FWHM for soft photons is wider than for hard photons; the 
smaller integration radius would make the spectrum harder, re- 
sulting in a systematically higher temperature. This is a serious 
effect for a distant cluster with apparent core size larger than 
1' when the GIS integration radius is smaller than ~ 4'. Note 
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that 6' corresponds to ~ 1 Mpc at z = 0.1, and it mostly covers 
the cluster region corresponding to A c = 500 (See ^5.2.1l for 
definition). The systematic error in estimating the bolometric 
luminosity due to the fixed integration radii will be discussed 
in the next subsection. 

We subtract background spectra that were obtained during 
blank-sky observations. Because the ASCA background has a 
detector position dependency, we extract them from the same 
region as the cluster in the detector coordinates. 

The instrument response can be split into two parts: a re- 
distribution matrix (RMF), which specifies the channel proba- 
bility distribution for a photon of given energy, and an effec- 
tive area curve (ARF), which specifies the telescope area and 
window absorption. We utilized the latest version of the GIS 
RMFs, gis2v4_0.rmf and gis3v4_0.rmf, while we generate the 
SIS RMF using the FTOOL sisrmg. We built the ARF files with 
the ASCAARF program appropriate for the cluster extended 
emission, summing the ARFs for each bin in the cluster image 
region according to the weight of the photon counts. 

4.2. Fitting with Raymond-Smith model 

In order to measure the average, emission-weighted X-ray tem- 
perature of the gas, we fitted the SIS and GIS spectra simul- 
tane ously with a thin-the rmal plasma emission model from 
Ravmond & Smith! Jl977l) . There are four parameters in the 
spectral model, the temperature kT, the metallicity relative to 
the solar abundance Z, the redshift z, and the normalization 
factor. The redshift of each object was fixed at the cataloged 
value in the NASA/IPAC Extragalactic Database (NED). The 
fitting function was convolved with the telescope and detector 
response functions. In the spectral fitti ng, we used the XSPEC 
version 9.0 analysis software (Arnaud, 1996). We rebinnedthe 
spectral channels so that each bin contains at least 40 photons. 

lYaaoobl d 19991) pointed out that fitting with fixed A^h re- 
sult in a systematically high temperature because of the seri- 
ous decrease of the low-energy efficiency of the SIS since early 
1994. To avoid the problem, the absorption column density A^h 
was allowed to vary. Note that in the observation of #57 A402, 
which was done during the A07 phase, a serious reduction in 
the efficiency below 1 keV was seen in the SIS spectra, so we 
used only the GIS data for this cluster. 

We checked the contribution of foreground/background 
contaminating sources in the ASCA spectra. We picked up the 
point sources whose photon counts are greater than 10% of 
the cluster from the HRI source lists. Excluding regions of 
r = 1' circles around the sources, we recalculate the spectra, 
responses and backgrounds to determine the temperature by 
the Raymond-Smith model fitting. Note that Energy Encircled 
Function at 1' is approximately ~ 0.3 for the GIS and ~ 0.4 for 
the SIS, respectively. We estimated the difference of the best-fit 
temperatures with and without point sources excluded relative 
to the lcr error of the measurement, and found that the contri- 
bution of the point-like sources does not affect the temperature 
measurement for most cases, except for #56 CL0500-24, #66 
CL0024+17, #76 3C220.1, and #20 A963. For A963, we ex- 
cluded one nearby source from the spectral region of the GIS 



and th e SIS. More spec ific analyses a r e desc ribed in lOta et all 
Jl998t) for CL05 00-24, ISoucail et alJ feOOOh for CL0024+17, 
lOta et alJ fcOOOh for 3C220.1. In the case of CL0500-24 and 
CL0024+17, only the SIS data were used. Table 4 lists the re- 
sults of the spectral analysis. 

We compared our results with values published in IWhitel 
( 2000) (there are 26 clusters in common) to find a good agree- 
ment between the two measurements with a mean temperature 
ratio of 1 .04. We also compared 6 high redshift c lusters with 
constrained ASCA temperature with Ettor i et al. ] i2004 . The 
temperatures for 5 clusters are statistically consistent between 
the two results, while there is systematic difference for #71 
3C295. As for 3C295. lEttori et al. I ( Eool excluded the cen- 
tral emission associated with the AGN, thus our temperature 
measurement may be affected by the AGN. 

We estimated the X-ray luminosity in the 2-10 keV band 
from the GIS flux, Lx(6'), and converted it into the bolomet- 
ric luminosity, Lx,boi(6'X using the emissivity of the Raymond- 
Smith plasma model. In order to check the systematic error of 
the luminosity estimation, we calculated the bolometric lumi- 
nosity by integrating the /3-model surface brightness distribu- 
tion determined with the ROSAT HRI within the same inte- 
gration area, to find that two estimations are consistent within 
about 15%. Then we derived the bolometric luminosity within 
r 5 oo (see section l5~2.1> . Lx,boi, by multiplying L x ,boi(6') with 
the ratio of the luminosities within rsoo to 6' using the j6-model. 
In Table 4 we listed Lx(6'), £x,boi(6'), and Lx.boi- 

5. Cluster parameters 

In ^3] we analyzed the cluster surface brightness distribution 
under the assumption that the gas is isothermal, and the spatial 
distribution is described by the single y6-model or the double 
/3 model. In |4j we determined the average X-ray temperature 
with the ASCA spectral data. The X-ray images and the fitting 
results of individual clusters are shown in Fig. IC.ll We will 
derive some physical quantities of the clusters from the X-ray 
parameters obtained from the analysis and investigate the red- 
shift dependence of these parameters and the histograms. 

5.1. Parameters from spectral analysis 

5.1 .1 . Redshift dependence of spectral parameters 

We show the spectral parameters obtained with the Raymond- 
Smith fittings as a function of redshift in Fig. ^3 in the plot 
of the tem perature , we quo t ed the results of nearby clusters 
( Mohr, Mathie sen & Evrardl Il999l) for comparison. There is 
no significant change in the temperature over a wide redshift 
range, particularly in z < 0.5. On the other hand, some clus- 
ters with z > 0.5 resulted in very high temperatures, though 
their errors are not well-constrained due to the limited photon 
statistics and the error range overlaps with the high temperature 
clusters in z < 0.5. 

As indicated from Fig. II Oh . it is hard to detect a clus- 
ter whose X-ray flux is below ~ 10~ 13 ergs^'cm" 2 due to 
the detection limit of ASCA. This flux corresponds to ~ 1 x 
10 44 erg s 1 at z = 0.5 and to kT ~ 2.5 keV from the L x - T 
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Fig. 10. X-ray temperature (a), metal abundance(b), the 2-10 keV flux(c) and the 2-10 keV luminosity within r < 6'(d) mea- 
sured with ASCA. At z > 0.1, the filled circles and the open triangles are the regular and the irregular clusters in our sam- 
ple. The filled stars show the double-/? clusters in our sample. The ASCA sensitivity curve is indicated with the dotted line 
in the panel of kT, assuming the luminosity-temperature relation. The temperatures of clusters with z < 0. 1 were taken from 
iMohr. Mathiesen & Evrard {l999). The asterisks and the filled squares denote the single-/? and double-/? clusters in the nearby 
sample, respectively. 



relation. Thus at z > 0.5, clusters with temperature lower than 
2.5 keV are barely detected. We consider that this can account 
for the higher average temperature at z > 0.5. 

5.1 .2. Histograms of the spectral parameters 

In Fig.^2 we plot the number of occurrence of each spectral 
parameter to study how the samples are distributed in the pa- 
rameter space, regardless of the redshift. We show the distri- 
butions of the regular and the irregular clusters separately as 
well as the distributions of all the clusters. As a result, there is 
no clear difference between the distributions of the regular and 
irregular clusters. In Table 8 we summarize the mean and the 
standard deviation of the spectral parameters. 

5.2. Parameters from fi-model analysis 

5.2.1 . Definitions of cluster parameters 

From the /3-model fits and spectral analysis we obtained four 
primary X-ray parameters (kT ,/3,r c ,S o) to describe the intra- 
cluster gas. From these parameters, we evaluate several impor- 
tant quantities that characterize properties of the clusters and 
the ICM. Below we summarize the definitions of these quanti- 
ties. 

1 . Electron density at the cluster center : n e o 

The central electron density can be determined from the 
central surface brightness So, /3, r c and the temperature. 
However, from the present analysis, we obtained the cen- 
tral surface photon flux convolved with the telescope and 
the detector responses. In this case So can be written as 

S p0 (£i,£2) = j dE' j dER(E',E)A{E)n M n m 
A p (T,Z,E,z)^r c T(3j3- 1/2) 



MD A (i+z)) 2 r(3/?) 



(3) 



where S p o(Ei,E2) [counts s^'crrT 2 ] is the photon flux 
in the E\ - E2 energy band, R(E',E) and A(E) 
respectively represent the response function of the 
detector and the effective area of the X-ray tele- 
scope/detector system, A p (T, Z, E, z) the emissivity in units 
of photons s _1 cm 3 keV _1 for the object at redshift z. Da 
is the angular size distance to the cluster. We utilized the 
XSPEC program to perform the convolutions with the de- 
tector response functions. We adopt «ho = (/W^H)«e0, 
where p.n = 1-40 and^ c = 1.167. 



2. Cluster mass profile and density profile: M(r), p(r), po 
Assuming that the intracluster gas is in hydrostatic equilib- 
rium, the following condition is satisfied: 



kT jdlnp„ as dlnT 



yomip \ d In r dlnr 



GM(r) 



(4) 



In the case that the cluster gas is isothermal and has a /?- 
model density profile, p gas (r) = p gas (0)(l + (r/r c ) 2 )~ 3 P /2 , 
where p gas (0) = p e m p n e o, the total mass contained within 
the radius r, is estimated from 



M(r) = 



3kT/3r (r/r c ) 2 



firripG 1 + (r/r c ) 2 ' 
The average density within r is then 
M(r) po 



P(r) 



(5) 



(6) 



where po = 9kTfi/4-7TG/j.m p r 2 is the central total matter den- 
sity. 

3. Cooling timescale of the gas at the cluster center : f coo i 
We estimate the radiative cooling timescale of the intraclus- 
ter gas at the cluster center as 



^cool — 



<7ff«e0 ' 



(7) 



where qg is related to the volume emissivity of thermal 
Bremsstrahlung through e$ = q$n 2 T l l 2 . 
Cluster limiting radius (overdensity radius) and cluster 
mass : ^500 and M500 

We determine a cluster limiting radius within which the av- 
erage density p(r) is equal to A t . times the critical density 
of the universe at the collapse time; namely 



p(r) = A e p crit (z co i). 



(8) 



We adopt a fixed overdensit y of A r _ 500, which is jus- 
tified in a sense that lEvrard. Metzler. & Navarrol |l996) 
suggested from their numerical simulations to use this 
value to study the gas properties and that the hydro- 
static assumpt ion is not valid beyond this radius, and that 
Finoauenov. Reiprich & Bohrineer J2001I) showed that an 
assumption of isothermality also works at such overdensi- 
ties. Since we do not know the redshift of the cluster col- 
lapse, the most simple assumption is that the clusters are 
observed just after they are formed, i.e. z co i = z b s . We will 
determine the overdensity radius, r^oo under this assump- 
tion and calculate the hydrostatic mass within rsoo, M500, 
from Equation|5] 
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Fig. 11. Histograms of the spectral parameters determined with ASCA for 79 clusters (open). The panels (a) and (b) show the 
best-fit temperature and the metal abundance determined from the Raymond-Smith model fitting, respectively. The X-ray flux 
and the luminosity within r < 6' in the 2-10 keV band estimated with the GIS are shown in (c) and (d), respectively. The results 
of the double-/? clusters are superposed (hatched). 

Table 8. Means and standard deviations of cluster parameters 



Parameter 


Regular (45) 


Irregular (34) 




All (79) 




Mean 


Standard deviation 


Mean 


Standard deviation 


Mean 


Standard deviation 


kT [keV] 


6.3 


2.8 


7.2 


2.7 


6.8 


2.8 


Z [solar] 


0.33 


0.22 


0.24 


0.15 


0.29 


0.19 


Lxjxa [erg/s] 


1.7 x 10 45 


1.9 x 10 45 


1.7 x 10 45 


1.7 x 10 45 


1.7 x 10 45 


1.8 x 10 45 




0.56 


0.11 


0.73 


0.61 


0.64 


0.32 


r c [h£ Mpc] 


0.076 


0.060 


0.273 


0.259 


0.163 


0.202 


rieo [hj'g cnT 3 ] 


3.6 x 10~ 2 


2.9 x 10~ 2 


7.4 x 10~ 3 


8.4 x 10~ 3 


2.4 x 10~ 2 


2.6 x 10~ 2 


Po [/i^gcm -3 ] 


2.5 x 10~ 24 


2.5 x 10~ 24 


7.3 x 10- 25 


2.7 x 10~ 24 


1.7 x 10~ 24 


2.7 x 10- 24 


fcod [Gyr] 


4.4 


4.6 


15.6 


8.1 


9.4 


8.5 


r 500 [Kjl Mpc] 


0.96 


0.22 


1.1 


0.4 


1.0 


0.3 


M 500 [/d Mq] 


4.3 x 10 14 


5.1 x 10 14 


9.7 x 10 14 


2.4 x 10 15 


6.7 x 10 14 


1.6 x 10 15 


M eas [h^ /2 M ] 


7.1 x 10 13 


4.3 x 10 13 


9.5 x 10 13 


6.5 x 10 13 


8.2 x 10 13 


5.5 x 10 13 


/ gas [^,o /2 ] 


0.20 


0.08 


0.18 


0.07 


0.20 


0.07 



5. Gas mass, and gas-mass fraction within rsoo : M gas and / gas 
The gas mass within rsoo is derived with 

J- 
p gas (r)4nr 2 dr 
o 

= 4wp gas (0)r c 3 (1 + x 2 y 3l3/2 x 2 dx, (9) 
Jo 

where x = r/r c and X500 = rsoo/rc- Then we obtain the 
gas-mass fraction with / gas = M g . ds /Msoo. 

Among the four X-ray parameters, the temperature kT is 
independently determined from the other three y6-model param- 
eters. However, the statistical errors of the three parameters are 
coupled to one another. In particular the coupling between j3 
and r c is strong (Fig. 0. We determined the statistical errors 
of the cluster parameters listed above with this coupling taken 
into account. For that purpose we first determined the error do- 
main, i.e. the statistically allowed parameter region, in the four 
dimensional parameter space. Then, evaluating the cluster pa- 
rameters for all combinations of the X-ray parameters in the do- 
main, we determined the maximum and the minimum param- 
eter values of the domain. For the double-/? clusters, we also 
calculated those cluster parameters from the double /3-model. 
The methods of calculation are similar to those shown above. 
The details are shown in AppendixlAl 

5.2.2. Redshift dependence of X-ray parameters 

We plot the /?-model parameters, and the parameters derived 
from those X-ray parameters as functions of redshift in Fig. 
IT21 We show the results of double-/? model fits, and their in- 
ner and outer components are distinguished by different sym- 
bols. In the figures we also plotted the parameters taken from 
Mohr. Mathies en & Evrard Jl999h for clusters with z < 0.1. 



We do not see in those figures any clear redshift dependence 
in the distributions of the X-ray parameters except for the pa- 
rameters related to rsoo shown in Fig. ^] We will go back to 
these parameters in section l5.2.3l and focus on the parameters 
that do not involve the overdensity radius. 

From Fig. 112b. we notice that the c ore radius shows anap- 
parent redshift dependence. As noted in IOta&Mitsudalll2002l) . 
the core radius shows a remarkably large cluster-to-cluster dis- 
persion, spanning over two orders of magnitude. The core radii 
of the irregular clusters are systematically larger than those of 
the regular clusters, and there seems to be a gap in the r c distri- 
bution at around 0.1 /i^Mpc. The regular clusters also show a 
similar bimodal distribution in r c but the fraction of the larger 
r c group decreases with increasing z. 

To investigate the selection effect, we created simulation 
clusters with the Monte-Carlo method and performed the anal- 
ysis on the simulation clusters. We scale the count rate of 
MS0906.5+1110 at z = 0.18 to estimate the expected total 
counts for a 40 ksec observation of a cluster at z = 0.5 with 
a typical luminosity Lj ~ 1 x 10 45 ergs -1 . The expected 
cluster counts are about 600 counts. We simulated a series 
of cluster images with various core sizes and found that the 
signal-to-noise ratio is quite low for clusters with large cores 
of r c > 400 kpc at z > 0.5. Based on the results we es- 
timated the sensitivity of the current HRI observation to be 
Sq ~ 3.6 x 10~ 3 counts s^'arcmkr 2 . We show the sensitivity 
curve in Fig. 112b . Thus the redshift dependencies of the core 
radius can be explained by a selection effect. 

Therefore we conclude that the X-ray parameters, tempera- 
ture, core radius, /?, and the central electron density are consis- 
tent with showing no significant trend of evolution at z < 0.5. 
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Fig. 12. Results from j6-model analysis. In the panels (a)-(e), /?, r c , n e , po, and f coo i are shown. At z > 0.1, the filled circles and 
the open triangles are the regular and the irregular clusters in our sample. The filled stars and the open stars show the inner core 
and the outer core of the double-/? clusters in our sample. At z < 0-1 the filled squares and t he open squares denote the inner 
and outer components of the double-/? clusters in the nearby Mohr. Mathiesen & Evrard ( 1999) sample. The asterisks denote the 
nearby single-/? clusters. In the panel of r c , the selection effects due to the sensitivity and the spatial resolution of the ROSAT HRI 
are indicated with the dotted line and the dashed line, respectively. In the panel of f coo i we show a curve on which f coo i is equal to 
the age of the Universe at the cluster redshift. 

Fig. 13. rsoo, M500, M gas , and / gas derived from /?-model analysis are shown in the panels (a)-(d). See 35.2.1l for definitions of the 
parameters. The meanings of the symbols are the same as in Fig. 1121 



5.2.3. Gas-mass fraction within r 50Q and the 
systematic error 

Although the X-ray parameters that are directly determined 
from observation do not show strong redshift-dependence, we 
find a weak redshift dependence in the overdensity radius, 
which approximately follows rsoo oc (1 + z o bs)~ ' 6 (Fig. I13t0. 
The dependence is likely to be introduced when we define 
rsoo by Equation |8j i.e. p(r) = A c p c (z bs)- In other words, the 
redshift dependency is introduced by the assumption, z co i = 
Zobs- We confirmed that such a negative dependence disappears 
when we assume a constant z co i independently of z b s , for ex- 
ample Zcoi = 1. 

In Fig.ll3H. we show the gas-mass fraction inside r^oo de- 
termined with z co i = Zobs- We do not see a clear dependence on 
the redshift in our own data. We then obtain the average gas- 
mass fraction of our sample of 79 distant clusters to be 

</ gas > = (0.20 + 0.07) h~l 12 , (10) 

where the quoted error is the standard deviation of the cluster- 
to-cluster variation. On the other hand the average gas-mass 
fractio n of the nearby samples of iMohr. Mathiesen & Evrard! 
(1999) is (/gas)nearby = 0.12 ± 0.03 (Note we recalculated 
the value with Equation E] and A r - 500 under t he ass ump- 
tion of Zcoi = Zobs)- IMohr. Mathiesen & Evrard ill 9991) esti- 
mated their systematic errors to be ~ 10%. Although our 
distant sample seems to show higher / gas values in compari- 
son to the nearby samples, the two results are consistent with 
each other within their errors. If we further compare our re- 
sult with the baryon density in the Unive rse determined by th e 
Wilkinson Microwave Anisotropy Probe ( Sperg el et all l2003). 
Qb/Qo = 0. 16, it is again higher although they are in agreement 
within the errors. We will thus examine the possible systematic 
errors of / gas estimation due to (1) the choice of limiting radius, 
(2) the calibrations of the X-ray telescope/detector systems, (3) 
the assumption of z co i = Zobs, and (4) the effect of the tempera- 
ture gradient. 

(1) We have defined the cluster limiting radius with A c = 
500 and the resultant value is typically r^oo ~ 1 Mpc for the 
current sample. On the other hand, the extent of the observed 
X-ray emission, rx, is found to be larger for most of the clusters 
(Fig. IC.lt and rx/r^oo = 1.5 on average. Thus we do not need 
to worry about the effect of extrapolation in the current / gas 
estimation. The average gas-mass fraction within rx is derived 
as (0.24 ± 0.06) /z~„ /2 , which agrees with Eq. [10| within the 
errors. 



(2) Since the determination of the gas-mass fraction re- 
quires absolute calibrations of the X-ray telescope/detector sys- 
tems of ROSAT and ASCA, we carefully examined the calibra- 
tional errors and found that they can cause at maximum 25% 
errors in / gas (see AppendixlElfor details). 

(3) The assumption of z co i = Zobs may be a source of un- 
certainty in / gas . We can infer z co i from the condition that the 
central mass density should be higher than the average mass 
density, namely po > p(rsoo)- Since the observed range of po is 
~ 2 x 10~ 26 - 1 x 10~ 23 g cm" 3 , we obtain z co i < 1-3 - 18, where 
the smaller (larger) value corresponds to the clusters with large 
(small) core radii. Thus it is likely that the clusters with large 
r c were formed at z co \ ^1-3. 

We then vary the formation redshift, z co \ for which we cal- 
culate the critical density, p C rit(z C oi) and find the radius where 
the measured matter density is 500 times p c rit(z C oi)- If we sim- 
ply assume a fixed formation redshift for all the clusters, rang- 
ing from 0.5 to 1.5, and calculate the mean gas-mass frac- 
tions within rsoo in the same manner as Eq^B (/gas) varies 
from 0.18 to 0.12 with a typical standard error of 0.05. Thus 
the / gas estimation largely depends on the assumption of z co i- 
However such an effect is expected to be more serious for 
the low-redshift clusters, whose a verage gas-mass fraction was 
measured to be <f oa s)nR a r hv ~ 0-12 (Mohr, Mathi esen & Ev rard. 
1999; Sander son et all 120031) . Thus we may not attribute the 
systematic error in measuring / gas for the distant sample to the 
assumption of z co i- 

(4) The emission-weighted temperature reflects the temper- 
ature of the cluster core region. Then, if there is a significant 
temperature drop at the center, it may cause an overestimation 
of the gas-mass fraction because the cluster hydrostatic mass 
estimation is more sensitive to the temperature profile than the 
gas mass. Such temperature drops were usually found in clus- 
ter cores with short (~ a few Gyr) cooling timescales. The spec- 
tral analysis of the cooling flow clusters with the XMM-Newton 
and the Chandra satellites showed that the temperature drops 
typically by a factor of 3 over the central r < 20 kpc re- 
gion (e.g. iTamura et all l200ll ISchmidt etafl l200lh . We then 
estimated the emission- weighted temperature within 1.5 Mpc 
(corresponding to the typical integration radius for the GIS 
spectra) assumi ng the radial temper ature profile of T(r) oc r° 2 
(from Fig. 1 of Ta mura et all ( 1200 II) ) and the /?-model surface 
brightness distribution with r c = 50 kpc and /3 = 2/3 to find 
that it is lower by about 30% than that of the outer(r > 0.1 
Mpc) region. On the other hand, we obtained the mean gas- 
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mass fraction for the 26 regular clusters with f coo i < 3 Gyr to 
be </ gas > = 0.22 ± 0.08, which is larger by 20% than that of 
the rest of the sample. Thereby we consider that, for at most 
1/3 of the samples, / gas may be overestimated to some extent, 
depending on the degree of the temperature gradient. However, 
considering the fact that the correlation between / gas and ? coo i 
is weak as well as that the ranges of / gas for clusters with short 
and long cooling timescales are not very different, the effect of 
the temperature gradient in estimating the mean gas-mass frac- 
tion for all the samples is suggested to be not large compared to 
the cluster-cluster variation. For more accurate measurements 
°f /gas, we need to constrain the temperature profiles for the 
individual clusters and reduce the measurement uncertainties. 

5.2.4. Histograms of the X-ray parameters 

Since we found no significant evolution in the X-ray param- 
eters, we will investigate the distribution of X-ray parameters 
obtained from the /3-model analysis without distinguishing the 
clusters by redshift. In Fig.^J we show the histograms of the 
X-ray parameters. We find that both kT and /? are distributed in 
ranges smaller than 1 order of magnitude, while r c and n e Q are 
distributed over almost two orders of magnitude. 

While the mean value of r c for all the clusters is (r c ) = 
0.163 /z^jMpc, if we treat the regular and the irregular clus- 
ters separately, we obtain (r c ) = 0.076 /i^jMpc and (r c ) = 
0.273 /i~ ( jMpc for the regular and the irregular clusters (see 
Table 8). Thus these are different by a factor of 3. We also 
notice that the distributions of the irregular and regular clus- 
ters are not separated. Instead, the distribution of the regular 
clusters has a double-peaked structure, whose core radii corre- 
sponding to the two peaks are 50 /z^'kpc and 200 /i^kpc re- 
spectively. Thus there is about a factor of four difference. The 
peak of the larger core radius coincides with that of the irreg- 
ular clusters. One may consider that the regular clusters with a 
large core radius were classified as regular because the count- 
ing statistics were not very good. However, we find that the 
statistics of irregular and regular clusters with a large core are 
not very different. Thus it is difficult to explain this coincidence 
just by statistics. The core r adius distribution when the nearby 
Mohr. Mathie sen & Evrardi lll999h samples a nd the our dis- 
tant sa mples were added together was shown in lOta & Mit suda 
i2002h . 

In Fig. I14b . we also show the distributions of the core 
radii of the double-/? clusters. We notice that the double-peaked 
structure of the core radius of the regular clusters seems re- 
lated to the double-/? structure because the core radius of the 
smaller-core component is distributed around the lower peak 
of the regular-cluster core radius distribution, while the larger- 
core is around the higher peak. 

In section l3~5l we have shown that about 20% of regular 
clusters have significant double-/? structures, and for about 60% 
of regular clusters, the existence of a similar structure cannot 
be rejected. We have also shown that there are inner-core dom- 
inant and outer-core dominant cases and that a single /?-model 
fit picks up the core radius of the dominant core (Fig.|9j. Thus 
we consider that the correspondence between the core radius 



distribution of all the regular clusters and the distribution of 
double-/? clusters is not just a coincidence but that it is related 
to the double core nature of the regular clusters. 

The electron density, the dark matter density and the cool- 
ing time at the cluster center show similar double-peaked dis- 
tributions. On the other hand, as shown in Fig.^] rsoo and the 
other three parameters evaluated within rsoo are distributed in 
ranges smaller than 1 order of magnitude. We will discuss the 
correlations between the parameters in the next section in more 
detail. We list the mean and the standard deviation for all the 
parameters in Table 8. 

6. Correlations between cluster parameters 

We investigate correlations of various cluster parameters with 
the X-ray temperature, kT in 36. H and also with the core radius, 
r c in 36.21 which showed the distinct double-peaked distribu- 
tion and thus may provide a clue to understand the structures of 
clusters. Possible systematic errors will be considered in each 
subsection. 

In the following analysis, we first calculate the 
Pearson's correlation coefficient, Rxy = Z/C^/ ~ _ 
?)/ VZ/(X; - X) 2 ^Zi(Yi - ?) 2 , to measure the strength of the 
correlation. We take the logarithm of the two measured param- 
eters, x and y, namely X = log x and Y = log y, except for the 
case of the Z - T and the fi—T relations. Then if \Rxy\ ^ 0.3, 
we derive the best-fit relation between the two parameters 
assuming the power-law function. In order to take into account 
the statistical uncertainties of both the x and y axes, we 
performed the x 1 minimization in the linear (Y = aX + b) fit 
by defining x 2 = Zi( Y i ~ ( aX i + b)) 2 /(acr 2 XJ + cr 2 .), where 
<Tx,i and cry, are the lcr errors for the parameter X, and T,, 
respectively. Because in all the cases below the fits were not 
statistically acceptable due to large scatters of the data at the 
90% confidence level, the error ranges of the coefficients a 
and b were estimated from the dispersions of the data points 
around the model functions rather than the photon statistics. 
We excluded two irregular clusters, #43 A1758 and #79 
MS1054.5-0321, in the analysis except for the Z - T and the 
Lx.boi - T relations because their /?-model parameters were not 
well constrained ( 33.4> . although they are plotted in Fig. [HO 
and ED We then compare some of the resulting relations with 
the predictions of the self-similar model. 

6.1. Correlations with the gas temperature 

InFig.H6ji~h, we show eight parameters derived from the spec- 
tral and the image analysis, the gas mass, and the gas-mass frac- 
tion etc. as a function of X-ray temperature. We mainly show 
the relations derived for A c = 500 below. In Table 7, we show 
the cluster parameters of the individual cl usters for the over- 
densi ties of A c = 500Q 427 and 18^ 2 fi 427 dNakamura & Sutol 
Il997l) and in Table 9 the scaling relations with and without con- 
sideration of the co smological factor, E, = H z /Hq as noted by 
lEttori et al~l J2004I) . Since we did not find any strong redshift 
evolution in the observed X-ray properties, we first derive the 
parameter correlations regardless of their redshifts. We will dis- 
cuss the correlations in the case where we take into account the 
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Fig. 14. Histograms of the X-ray parameters determined with the single /?-model for 79 clusters (open). The results of the double 
/?-model fittings for the nine double-/? clusters are superposed in the panels, where the hatched and filled regions show contribu- 
tions of the inner core and the outer core, respectively. 

Fig. 15. Histograms of rsoo, M500. M gas and / gas derived with the single /? model for 79 clusters (open). The results of double /? 
clusters are superposed (hatched). 

Fig. 16. Relations of the metallicity (a), the bolometric luminosity (b), /? (c), the central electron density (d), the overdensity 
radius (e), the cluster mass (f), the gas mass (g), and the gas-mass fraction (h) with the X-ray temperature. In all the panels the 
results of the single /?-model fitting are plotted and the single-/? regular clusters, the single-/? irregular clusters, and the double-/? 
clusters are denoted with the filled circles, the open triangles, and the filled stars respectively. The horizontal and the vertical 
error bars are lcr. In the panels (b) and (e)-(h), the best-fit power-laws for the entire sample are shown with the solid lines. In 
the panels (a), (c) and (h) the sample means are indicated with the dashed lines. In panel (b), the best-fit L - T relations for two 
subgroups with r c < 0. 1 Mpc and r c > 0. 1 Mpc are indicated with the dashed and dash-dotted lines, respectively (see section l6~Tl 
for details). In panels (e)-(g), the slopes expected from the self-similar model are shown with the dotted lines. 



sample redshifts later. Note that kT and r c are in units of [keV] 
and [fe^Mpc], respectively. 

-Z-T 

In Fig. I16f a). there seems to be a slight decline of the 
metallicity, Z, against the temperature, however, we found 
the correlation coefficient, Rxy to be -0.23 for 78 clusters, 
where #66 CL0024+17 was not used in the fit because Z 
was so uncertain t hat the value was fix ed at 0.3 solar in 
the spectral fitting ( Sou cail et all I20001) . Thus there is no 
clear Z-T correlation. This is consistent with the previ- 
ous results on the metallicity measurements fo r 21 clusters 
reported by Mushotzkv & Loewenstein ( 1997). We thus in- 
dicate the mean metallicity of our sample, (Z) = 0.29, with 
the dotted line in the figure (see also Table 8). 
- Lx,boi - T 

We find a significant correlation between the bolometric lu- 
minosity and the temperature, and Rxy is 0.59 for 77 sam- 
ples. Note that #56 CL0500-24 and #66 CL0024+17 are 
not included in order to avoid the systematic difference in 
the integration radius for the spectra because the only SIS 
spectral data were used for these two clusters due to the se- 
rious contamination from the point sources in the vicinity 
of the clusters (see |4j. We thus obtain the temperature- 
luminosity relation for the 77 clusters to be 

Lx,boi [ergs- 1 ] = 6.53^ x 10 42 (£r) 3 08 --. (11) 

X 2 /d.o.f = 1703/75. As a result, the slope of the relation 
is consistent with those previous l y published for the nearb y 
clusters (e.g. iMarkevitchl U998t lArnaud & EvrardI fl999). 
However, if compared to the relation for higher redshift 
clusters bv lEttorietaLlfe004 . they showed a steeper slope 
of 3.72 ± 0.47 (lcr error). 

The large error in the normalization factor in Equation II II 
is affected by the intrinsic scatter around the best-fit rela- 
tion. We found from Fig. 116b that the normalization factors 
of the Lx.boi - T relation are significantly different between 
the regular and the irregular clusters: for a fixed tempera- 
ture, the regular clusters tend to have a larger luminosity 
in comparison to the irregular clusters. The difference be- 



comes more evident if we divide the sample into two sub- 
groups according to the core radius ranges of r r < 0.1 hzl 
Mpc and r c > 0.1 h~l Mpc (see also lOta & Mitsudal 12002). 
If we fit them separately, fixing the slope at the best-fit value 
of 3.08, we obtain 

iotboi [ergs" 1 ] = 9.12!^ x 10 42 (£r) 308 forr c < 0.1,(12) 
L x ,boi [ergs- 1 ] = 3.98!°; 22 x 10 42 (kT) 3M forr c > 0.1,(13) 

and * 2 /d.o.f = 810.9/38 and 199.4/37, respectively. Thus 
the the normalization factor for the small core clusters is 
larger than that for the large core clusters at the 10cr signif- 
icance level. 

- p-T 

There seems to be a trend of larger /? values fo r higher 
temperatures, as previously noted by Schindler (1999). 
However Rxy = 0.27 for 77 clusters, hence the correlation 
is not clear. 

- na-T 

We see from Fig. I16H that there is a significant scatter of 
the central electron density, n e o. We obtained Rxy - -0.28 
and thus the correlation is not clear from the data. On the 
other hand, we notice n e o is strongly correlated with the 
parameter r c , which will be discussed in the next subsection 
in more detail. 

- >*500 - T 

A strong correlation between the overdensity radius and the 
temperature is found {Rxy = 0.82), whose best-fit relation 
is: 

r SO0 [h^Mp C ]=0.38^(kT) - 5 ^, (14) 

and* 2 /d.o.f is 151.2/75. The power-law slope of 0.53±0.04 
is consistent with the value predicted from the self-similar 
model, 0.5, within the error range. 

- M 50 o - T 

There is a tight correlation between the cluster mass and 
the temperature {Rxy = 0.88). The best-fit relation is 

M 50 o [h^M G ] = I.64-;; ;; x 10 13 (£r) L68 -S:n, (15) 

and* 2 /d.o.f is 206.5/75. The slope of 1.68!°}° is slightly 
steeper than that expected from the self-similar relation, i.e. 
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£1 M = 0.3, Q A = 0.7, Zcoi = z b s 




Relation 


A c = 500 




a c = sooaP*" a c = i^a^ 1 

o„o+U.4« . „„., -,n, +11.47 





L x ,boi - T 6.53!°-^ x Vf i (kT) )m -Mi 7.94" f„ x 10 42 («y ul -<w2 1.07+' '° x W 4i (kTY 



M-T 1.64+°^ x io 13 (ytr) L68 -a!? 2.07!°-| x lO 13 ^) 1 ' 68 ^! 3.31^g x l0 n (kT) L69 ^ 

M im - T 2.88+J^ x 10 12 OfcD L85 -^ 4.32« : g x 10 12 (/tr) L80 -a!^ l.l0 + ^» x lO'W) 1 - 67 -^ 

/ gas -r o.28;^(^r)- - 08 -a!4 o.33^(«y -' l! ° : " o.47^(^r°-' 9 -"° 

£:'L xio , - r 5.50^ : « 2 x io 42 (fcr) 3 - 12 -li:3s 7.24^ x io 42 (/tr) 3 - 00 -o:45 1.01^ x io 43 (/tr) 2 - 8<, -o:S 

E Z M- r 1.73^ x 10 13 (yfcr) L71 -a!i 2.11^ : « x 10 13 (/tT) L73 -a!! 3.51 + °g x 10" (W) 1 ' 71 -^ 

£ ; M g;ls - T 3.05;^ 4 x 10 12 (£r) L88 -^ 4.62+; ^ x lO 12 ^) 1 - 82 -^ 1.19^1^ x lO'W)'- 68 -^ 



M500 °c r . However, from Fig. I16F . we consider that 
the degree of the departure from the self-similar relation 
is marginal if compared to the typical size of the statistical 
error bars of the data points. 

We compare the M500 - T relation with those de- 
rived for two samples of nearby cluste rs (Table 1 of 
Fino guenov. Reiprich & Bohringer (2001)). Because they 
calculated it under a different cosmology: £2o = 1 and 
Hq = 50 kms^'Mpc -1 , we calculate the relation using the 
same set of cosmological parameters and obtain M500 = 
2.04+°' 11 x l0 l3 (kT) h65 -°™ for the distant clusters. In com- 
parison to their flux-limited sample (HIFLUGCS) and the 
sample with temperature profiles, the slope is in a good 
agreement within the errors but it should be noted that the 
best-fit normalization is about 30% smaller for the distant 
sample. This may be attributed to the significant redshift 
dependency of the critical density and will be worth further 
investigation in the light of the cluster formation redshift. 
We also showed the E Z M - T relation calculated for A c = 
500Q 427 in Table 9. The result is within a range consistent 
with the relation for z > 0.4 (Etto rietal. 112004 . 
M gas - T 

The gas mass within rson, M gas , is strongly correlated to the 
temperature (R X y = 0.61). The fitting gives 

M s , s [/^ /2 M ] = 2.88!°;^ x lO 12 ^) 185 ^, (16) 

and^- 2 /d.o.f is 548.1/75. Taking into account the current 
statistical errors, we found a marginal steepening of the re- 
lation in comparison to that predicted from the self-similar 
model, M gas K T 15 (see Fig. 116b ). The measured slope is 
also found to be consistent with that of the M500 - T relation 

derived above within t heir errors. 

For nearby clusters, iMohr. Mathiesen & Evrard ( 1999) 
obtained the relation to be M gas = (1.49 ± 0.09) X 
10 14 (£776 keV) L98±018 under Q = 1 and H = 
50 kms^'Mpc -1 , which is significantly steeper than the 
theoretically expected slope of 1 .5. On the other hand, 
IVikhlinin. Forman & Jonesl dl999l) reported based on the 
ROSAT PSPC data analysis of nearby regular clusters, a 
flatter relation in the form of M gas oc T l • 71±013 j utilizing 
a different method in determining the limiting radius (they 
defined the baryon overdensity radius of /?iooo> correspond- 
ing to the dark matter overdensity of ~ 500). 



If we calculate the distant M gas - T relation within 
rsoo using the same set of cosmological parameters, 
M gas = 4.57+j-^ x 10 12 (fcr)' 8 6 -°^. The relation is slightly 
flatter than that derived by IMohr. Mathiesen & Evrardl 
Jl999l) but within a range consistent w it h the re- 
sult of either IMohr. Mathiesen & Evrardl i 19991) or 
IVikhlinin. Forman & Jonesl i 19991) under the current 
measurement errors. 

As shown above, although the relation within rsoo obtained 
for our distant sample is consistent with the local relations, 
it is found to be le ss steep than that f ound for z > 0.4, 
E z M m oc t 237±0A7 ( lEttori etal.ll2004 . which is not con- 
flict with the view that the lower redshift clusters contain 
more gas for a fixed temperature, as pointed bv lEttori et al. I 

— /gas — T 

We obtained a small correlation coefficient, Rxy = -0.34. 
The power-law fitting resulted in 

/ g as[A v - 3/2 ] = o.28!^^r)- -<:K, (17) 

and;^ 2 /d.o.f = 349.8/75. Because the resultant slope in- 
cludes within the error, we conclude that there is not a 
significant temperature dependence in the distant sample. 
Thus we show the mean gas-mass fraction of the sample, 
</ gas > = 0.20,inFig.[T6h. 



In the above analysis, we have not included the effect of the 
cluster redshift. However, because a weak redshift dependence 
is seen in the overdensity radius, which is rsoo oc (1 + z O bs) _0 ' 6 
(see 35.2.3I I. we checked how the parameter correlations in 
equations 11414171 will be changed if we divide the sample into 
the low-z (0.1 < z < 0.3) and the high-z (0.3 < z < 0.82) 
subsamples. As a result, there are no significant changes in 
the scaling relations compared to equations 1 1 4H 1 71 except that 
the M gas - T for the high-z subsample resulted in a marginally 
steeper slope of M gas = 5.4;^ 2 x 10 11 (kT) 2 - 61 -™* Cf 2 /d.o.f = 
170.2/30). In order to place firmer constraints on the scaling 
relations, we suggest that it is important to gather more obser- 
vational data with higher sensitivities and also reexamine the 
assumption of the isothermal gas distribution and the forma- 
tion redshift, z co i as already pointed out in 35.2.31 
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Fig. 18. Relations of /? (a), the central electron density (b), the cooling timescale (c), the temperature(d), and the overdensity 
radius (e) with the core radius. The meanings of the symbols are the same as Fig.^] The error bars are l<x. In the panels (a)-(e) 
the best-fit power-laws obtained for 77 distant clusters are shown with the solid lines. In the panels (b) and (c), the best-fit power- 
laws for log r c < - 1 and log r c > — 1 are also shown with the dashed lines and the dot-dash lines, respectively. In the panel (e) 
four dotted lines correspond to four different constant values of rsoo/r c . 



Fig. 17. Venn diagrams which illustrate the relation between 
the X-ray morphology and the optical morphology (a) and the 
X-ray core radius and the optical morphology (b) (see sec- 
tion |^] for definition of the cD cluster.) In the panel (b), we 
divided the sample into three subgroups: small core (r c < 0.1 
Mpc) single-/?, double-/?, and large core (r c > 0.1 Mpc) single- 
/? clusters. 

6.2. Correlations with the core radius 

- (Optical morphology)-r f 

First we investigate the relation to the optical morphology 
of the clusters. Some clusters contain a central dominant 
elliptical galaxy, i.e. a cD galaxy. We refer to such clusters 
classified as Bautz-Morgan types I and I— II as "cD clus- 
ters". We looked up the BM types of our sample clusters 
in the NED database and showed the relation between X- 
ray morphology and the BM type, and the core radius and 
the BM type in Fig. I17h and 1 17b respectively. In Fig. 117b 
we divided the clusters into three subgroups: the small core 
single-/?, the large core single-/? and the double-/? clusters. 
We find that all the cD clusters are regular clusters and the 
clusters with a small core (i.e. small core single-/? + double- 
ts) tend to contain a cD galaxy. However, not all the small 
core clusters have cD galaxies, and the fraction having a 
cD galaxy is 36% (5 of 14). Therefore it is not a simple 
one-to-one correspondence. In nearby clusters, the typical 
X-ra y core radius of c D galaxy is measured to be ~ 10 
kpc (Ike be et all II 999) and is significantly smaller than 50 
kpc. Thus it is unlikely that the small core represents the 
potential distribution of the cD galaxy itself, though some 
connection may be possible. At present, however, the data 
of the optical morphology is available for only 47% of the 
sample. Thus in order to clarify the correlation between the 
central galaxy and the formation of the small core, we need 
to collect more optical data. Though the above discussion 
was based on the BM-type classification, we suggest that 
it is also meaningful to take into account the existence of 
giant ellipticals at the cluster center. 

- P - r c 

There is a weak trend of larger/? for larger r c (Rxy = 0.62). 
The fit yields 

£ = 0.73^r° H -° : °\ (18) 

and^f 2 /d.o.f is 711.1/75. However the fit is not statistically 
acceptable due to the huge x 1 value. Regarding this, we 
have to be careful about the parameter coupling in the /? 
model fitting because the correlation seems to follow the 
direction of the coupling (Fig.|5Jl- This occurs noticeably at 
r c > 0.1 Mpc. The correlation is strong for the large core 



clusters with r c > 0. 1 Mpc (Rxy = 0.72) while it is weak 
(Rxy = 0.09) for the small core clusters with r c < 0.1 Mpc. 
In addition, all the clusters that exhibit extremely large core 
radii (r c > 0.4 Mpc) are irregular clusters and their sur- 
face brightness distributions are highly inhomogeneous or 
bimodal (Fig. IC.H . Accordingly the current spherical /?- 
model can cause the tight /? - r c coupling particularly for 
the irregular systems. 
- n e0 - r c 

As is clear from Fig. 114b. the central electron density, n e o, 
has a double-peaked distribution, similar to r c . We see from 
Fig. 118b that there is a strong correlation between n e o and 
r c . The correlation coefficient is Rxy - -0.85 for 77 clus- 
ters. From the^ 2 fitting, we obtain 

n e0 [h^cm- 3 ] = 0.89!°; 3 ° x l0- 3 r~ 139 ^". (19) 

^ 2 /d.o.f is 1 185/75. We also find that the slope tends to be 
steeper for the small r c clusters: if we fit the data points for 
two different r c ranges separately, we obtain 

n e0 [/4o 2 cirT 3 ] = 0.13!^ x lO^r' 1 '* 1 ^ for r c < 0.1(20) 

n e0 [h^crxT 3 ] = 1.32+^° x KTVJ 1,10 ^ for r e > 0.1(21) 

;f 2 /d.o.f are 650.3/39 and 373.4/34, respectively. Thus the 
gas distribution in clusters with r c < 0.1 Mpc is concen- 
trated more than expected from the relation for clusters with 
r c > 0.1 Mpc. This may suggest that the small core and 
the large core components have different physical natures 
and/or they are at different stages of evolution. 

'cool 

According to Figs. 112b andll8b. f coo i is significantly shorter 
than the age of the Universe for the small core clusters and 
then the radiative cooling is suggested to be important. We 
find a very tight correlation between f coo i and r c (Rxy - 
0.87 for 77 clusters). The f coo i - r c relations are derived to 
be 

fcooi [yr] = 9.55!|^ x 10 10 r' 3l! °°' for 77 clusters, (22) 

fcooi [yr] = 31.6!^;i x lO 10 ^ 68 ^' for r c < 0.1, (23) 

fcooi [yr] = 3.63^; 2 | x 10 10 rf° !§:!l for r c > 0.1. (24) 

^/d.o.f are 690.0/75, 402.4/39 and 175.9/34, respectively. 
Thus one possible interpretation of the small core may be 
that the small core radius does not reflect the shape of the 
gravitational potential but that it reflects the cooling ra- 
dius inside which the X-ray emission is enhanced. However 
we consider this is unlikely for the following reasons. If 
the small core radius reflects the cooling radius, it should 
evolve with time. However we do not find a strong redshift 
dependence in the core radius. Moreover we estimated the 
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cooling radius at which the cooling time of the gas becomes 
equal to f age to find it is larger than 50 kpc for most of the 
small core clusters. 

- T-r c 

Furthermore, as shown in Fig. I18H . we find that the tem- 
perature does not show a clear core radius dependence. We 
obtained a moderate correlation coefficient, Rxy = 0.45. 
The best-fit relation is derived as 

kT = 13.7!|fr" - 007 , (25) 

however, the power-law fit was quite poor (^ 2 /d.o.f = 
2228/75). Since the emission-weighted temperature re- 
flects the temperature of the cluster core region, the ap- 
parent lack of any strong T — r c correlation suggests that 
the temperature gradient in the cooling region is not very 
large. This is consistent with the results from the XMM- 
Newton observations, which revealed that the temperature 
gradient of nearby cooling-flow clusters is smaller than 
that expected from t he standard cooling fl ow model (e.g. 
Tam ura et alll200ll) . |Peterson et alJ J 2003) noted that there 
is a significant deficit of emission with temperature lower 
than Tq/3 {To is the ambient temperature) in the RGS spec- 
tra. Given that the temperature profile obeys T(r) oc r° 2 , 
the second term in Equation|4]is negligible compared to the 
first term at r > 0.3r c for /? = 2/3. Thus except for the 
central r < 0.3r c region, the cluster mass profile can be ap- 
proximated with the assumption of a constant temperature, 
suggesting the gas density profile obtained from the current 
isothermal /?-model analysis reflects the underlying clus- 
ter potential distribution. We further investigate the relation 
between rsoo and the core radius below. 

- rsoo - r c 

Since the overdensity radius is determined almost indepen- 
dently from r c and approximately rsoo 00 T 1/2 , the above 
results indicate that correlation between rsoo and r c is much 
weaker than the expectations of the self-similar model. It is 
clear from Fig. I18fe that the distribution of the data points 
are inconsistent with the curves of r5oo/r c = constant. The 
observed range of r5oo/r c is about 2—40. The correlation co- 
efficient is Rxy = 0.58 for 77 clusters. Then the x 1 fitting 
gives 

r 50 o [h^Mpc] = 1.51^;^r°' 15! ° S for 77 clusters, (26) 

and^ 2 /d.o.f. = 954.5/75. We notice that the observed dis- 
tribution on the rsoo _ r c plane is much flatter than the 
curves for constant rsoo/r c values or it is rather concen- 
trated around the two peak values of r c . The departure from 
the self-similar relation is more prominent for small core 
clusters with r c < 0.1 Mpc. The correlation coefficients are 
Rxy = 0.21 and 0.36 for r c < 0.1 and r c > 0.1, respectively. 
Thus for the large core clusters, the best-fit power-law rela- 
tion is 

rsoo [*^Mpc] = 2M + ^f/ c 31 -°" for r c > 0.1, (27) 

and ^ 2 /d.o.f. is 108.0/34. If we further restrict r c to a 
very narrow range of 0.1-0.2 Mpc, we find a steeper slope 

7 1 +0.30 

of r 500 = 4.03^; 22 r c ' -° 26 (y 2 /d.o.f = 21.2/18). Thus 



those twenty clusters may satisfy the self-similar condition, 
r^oo oc r c . However, we suggest from eauations l26H27l that 
it is difficult to explain the formation of the cores, particu- 
larly for the small core clusters, by the standard picture of 
the self-similar model. 

6.3. Implications on the origin of two core scales 

From the above discussion, it seems difficult to explain the 
small core size that we discovered in the histogram either by 
the potential structure of the cD galaxy or the cooling radius. 
As long as we rely on the hydrostatic assumption and the ft 
model, the X-ray surface brightness distributions are likely to 
represent the gravitational potential structures of the clusters. If 
this is the case, the double-/? nature of the X-ray emission pro- 
file reflects the shape of the gravitational potential of the dark 
matter, which is likely to have two preferable scales of ~ 50 
kpc and ~ 200 kpc. 

In an effort to constrain the physical status of ICM in 
the dark matter potential, comparing the high-resolution X- 
ray observations to gravitational lensi ng observ ations will pro- 
vide another powerful test (e.g. lHattori et all Il997l) . Thanks 
to the improvement of spatial resolution achieved by Chandra, 
now measurements on the cluster mass profile down to < 5 
kpc scale are possible at such high redshifts (z ~ 0.3) and 
several authors have measured the dark matter distribution 



in the lensing clusters under the hydrostatic hypothesis (e.g. 


Arabadiis et al., 


2002; Xue & Wu , 2002; Ota et al., 2004). For 


example, lota et all 


(2 


004) showed from the high-resolution 



Chandra data of CL0024+17 (z = .395) and the com pari- 
son with the detailed lens modeling bv lTvson et alJ (fl998) that 
the cluster density profile is well reproduced by the double-/? 
model and the inner core also reflects the underlying dark mat- 
ter potential. They also noted that the core structure may be 
related to the past merging event as inf erred from the optical 
observations JCzoske et allkOOliEooll) . 

Recently Havakaw a et alJ J200 4) estimated the dark matter 
distribution in a nearby non-cD, regular cluster, Abell 1060, 
from the Chandra data analysis, without explicitly using the 
double yS-model, and found a central mass concentration at 
r < 50 kpc. Their result also supports the idea that dark matter 
may preferentially be accumulate d within a radius of ~ 50 kpc. 
On the other hand. lEttori et al. 1 12004) suggested that no sig- 
nificant double structure is seen in the high redshift sample. 
Thus considering from the above, the double-/? nature of the 
ICM discovered in the present sample may be much related to 
the history of the past merging and the relaxation process. The 
small core component might be attributed to the presence of 
dark matter subhalos due to the cluster merg ers or the internal 
structures in clusters (e.g. Fuii ta et all 12002). 

Furthermore, since there is clearly a tight coupling between 
the core radius and the radiative cooling time as shown in sec- 
tion the detailed treatment of the thermal evolution of the 
ICM will also be important. A number of numerical simula- 
tions including non-gravitational effects such as radiative cool- 
ing and galaxy feedback have been carried out and thus pro- 
vide a clue to the underlying physics in the cluster core regions. 
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However, the authors pointed out difficulties in regulating the 
central over-cooling and produce a constant-de nsity core (e.g. 
Pearceetal .1 1200(1 Im asai & Kitavamal J2004h recently pro- 
posed a quasi-hydrostatic model, which predicts a characteris- 
tic temperature profile with an asymptotic temperature for the 
central region being ~ 1/3 of the non-cooling outer region, as 
observed in nearby "cooling flow" clusters. Thus detailed com- 
parison of the X-ray data with their model regarding the tem- 
perature and density profiles will be important to understand 
the evolution of the ICM structure. We need further investiga- 
tions to put a stronger constraint on the origin of the double 
nature of the cluster structures, which is however beyond the 
scope of the present paper and will be discussed in a separate 
paper. 

7. Summary 

We have analyzed the ROSATHRl and the ASCA GIS/SIS data 
of 79 clusters of galaxies at redshifts of 0. 1 - 0.82 in a uniform 
manner. We determined the X-ray surface brightness profile 
from the ROSATHRl data utilizing the /?-model and the average 
temperature and the luminosity from the ASCA data. We found 
that the clusters can be divided into two subgroups, regular and 
irregular clusters, from analysis to determine the centroid posi- 
tion of the X-ray image. We then performed a statistical study 
of the X-ray parameters and investigated the trends for redshift 
evolution and the scaling relations against temperature and core 
radius. The major results are summarized as follows. 

1 . We did not find significant redshift evolution in the X-ray 
parameters of clusters compared to the nearby clusters: the 
temperature kT, the core radius r c , /?, and the central elec- 
tron density n e o at z < 0.5. 

2. Among the X-ray parameters, the core radius shows the 
largest cluster-to-cluster variation. The core-radius distri- 
bution shows two distinct peaks at 50 kpc and 200 kpc. 
For 20 % of the regular clusters, inclusion of a second 
/?-model component significantly improved the x 1 values 
of the surface-brightness fitting. We find that the two core 
radii of the double yS-model are distributed in relatively nar- 
row ranges consistent with the two peaks of the single-/? 
clusters. There is no significant evolution in the fraction of 
double-/? clusters within the observed redshift range. 

3. We investigated the correlations between the temperature 
and the cluster parameters including the spectral and the /?- 
model parameters, the cluster mass, the gas mass and the 
gas-mass fraction etc. For the M500 - T relation, we found 
that the power-law slope of 1.68^'jj is marginally steeper 
than that expected from the self-similar model but in a good 
agreement with the results for the nearby clusters, while 
the normalization factor is about 30% smaller for the cur- 
rent distant sample compared to the nearby sample. We ob- 
tained the M gas - T relation to be M gas oc (kT) l - S6 ^ for the 
overdensity of A t . = 500 and found a marginal steepening 
of the relation in comparison to the self-similar model un- 
der the current statistics. The / gas - T relation is found to be 
consistent with having no correlation with the temperature. 



4. We studied the parameter correlations against the core ra- 
dius. We found that only 36% of the small core single-/? 
and the double-/? clusters are cD clusters and thus it seems 
difficult to explain the presence of the small core by the 
cD potential itself although there may be some causal link. 
There are strong n& - r c and f coo i - r c correlations and the 
slopes tend to become steeper for r c < 0.1 Mpc. On the 
other hand the fact that there is not a clear T - r c correla- 
tion suggests that the temperature gradient is not large even 
in clusters with short cooling timescales, which is consis- 
tent with the XMM-Newton and Chandra observations of 
the nearby clusters. Thus as long as we rely on the hydro- 
static condition and the /3-model, our result indicates that 
the dark matter distribution is likely to show two preferable 
scales of 50 kpc and 200 kpc. 

5. We showed that the r^oo - r c relation derived from the X- 
ray analysis does not agree with the expectations of the 
self-similar model, suggesting that the assumption of self- 
similarity is not valid in describing the density profile of the 
ICM, particularly for clusters with small core radius. 

6. We obtained the average gas-mass fraction within 7-500 to 
be </ gas > = (0.20 ± 0.07) hZ^ 1 for the distant sample. The 
calibrational error is estimated to be about 25%. The cur- 
rent estimation is based on some simplified assumptions 
(for example, the isothermal gas distribution), which will 
be refined in future studies. 
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Appendix A: Mass profile of double-/? clusters 

Suppose that the density distribution of intracluster gas is char- 
acterized by superposition of the two jfi-model gas profiles. 



Pgas 



(r) = J^PgasM 



i=\ 



1 + 1^ 

n 



-3A72 



(A.l) 



Then the X-ray surface brightness distribution is given by inte- 
grating the X-ray emissivity along the line of sight, 



5( ^/^~/^ n,(r)2A(r ' Z)) ^' 

2 



1 + 



-3A+1/2 



(A.2) 



We refer to the second line of the above equation as the dou- 
ble /J-model and utilized in 33.51 We also obtain the thermal 
pressure of the gas: 



P{f)=Y u n i {r)kT. 



(A.3) 



(A.4) 



Then the total cluster mass is estimated from the hydrostatic 
equation to be 

kTr 2 d\nn(r) 

M(r) = - — , 

ymi p (j or 

£ al 1 b.html ^Z ; A»o,rr 2 [l + (^) 2 r 3ft/2 - 1 

pm p G n i( r ) 

The average density profile of the total cluster mass is then 

M(r) 9kT E,A»o,,rr 2 [l+(i) 2 ]- 3 A/ 2 -> 

= 4 7 = ~ A F v — rr ■ (A.5) 

±nr 3 4-npm p G 2_; n,-(r) 

We define the central cluster mass density as 

9kT IsiPinojrJ 2 



P0 = P(r)\r^0 = 



4-nfmipG n 0,i 



(A.6) 
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Appendix B: Systematic error of gas-mass fraction 

The gas-mass fraction is determined from the parameters, T, 
P, r c , and n e o, and n e o is further determined from the central 
surface brightness S p o, T, r c , and /3. Among those parameters, 
a possible systematic effect for r c was investigated in section 
13.41 and was found to be less than the pixel size of the X-ray 
image, ft is strongly coupled to r c and its systematic error is de- 
termined by the systematic error of r c . On the other hand, n e a 
requires the absolute calibration of the X-ray telescope/detector 
effective area, which usually contains large systematic errors. 
The temperature, T, requires a calibration of the effective area 
as a function of X-ray energy. Thus we will investigate the sys- 
tematic errors in n e Q and T below. 

Systematic errors in the electron density mainly come from 
the calibration of the ROSATHS1. Because the gain of the HRI 
continuously decreased from the launch to the end of the mis- 
sion, the conversion factor from the photon energy to pulse 
height is time dependent. Thus the errors in the gain deter- 
mination cause significant error in the absolute flux. The pos- 
sible range of th e gain variations over t he lifetime of ROSAT 
was measured by Prestwitch et al . (1998). Thus we tried the re- 
sponse matrices for the two extreme gain values in calculating 
Equation[3]to estimate the maximum systematic errors. Taking 
into account that the result also depends on the cluster emission 
spectrum, the effect is ~ 20% in the worst case. The calibration 
of the effective area of the X-ray telescope and the HRI system 
is also reported by the ROSAT Science Data Center, from which 
we estimate that the systematic errors in n e o is ~ 10%. 

Systematic errors in the temperature come from the calibra- 
tion of ASCA. The response functions of the ASCA XRT/GIS 
and XRT/SIS are well calibrated for point sources. However, 
for the extended sources, there still are significant systematic 
errors. In particular we found that the derived temperature is 
dependent on the spectrum integration region on the detector 
if the integration radius is too small. In the spectral analysis, 
we determined the radius so that the dependence becomes in- 
significant. However, we still consider there is some systematic 
effect related to this problem and estimate that it will affect / gas 
by about 5% (notice T affects both M500 and M gas ). 

In total, the systematic error due to the instrument calibra- 
tions is estimated to be ~ 25%. 



Fig. CI. ROSAT HRI images, radial surface brightness profiles 
and the ASCA SIS and GIS spectra of 79 clusters. The HRI im- 
ages are smoothed by a Gaussian filter with cr = 7". 5 - 15" and 
the contours whose levels correspond to n(= 3, 5, 9, 15, 31, 63) 
times the Icr background level are overlaid. The backgrounds 
are not subtracted. The X-ray centroids that were determined 
and used to derive the radial profiles in the image analysis are 
marked with the crosses. In the central panels, the crosses de- 
note the observed radial profile of the HRI, and the step func- 
tions show the best-fit fi models. The best-fit background levels 
are shown with the dashed lines. For the nine double /? clus- 
ters, the results of the double (3 model fitting are shown instead 
of the single (3 model and the inner and the outer components 
are also shown with the dash-dot and the dotted lines, respec- 
tively. The X-ray significance radius, r A , and the overdensity 
radius, rsoo are shown with the vertical dashed, and dotted lines 
respectively. In the right panels, the ASCA spectra fitted with 
the Raymond-Smith model are shown, where the crosses and 
crosses with circles denote the spectra obtained with the GIS 
and the SIS respectively, and the stepped lines show the best-fit 
models convolved with the telescope and detector responses. 
The fitting residuals are also shown in the panels. 
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